Solution of two— center time— dependent Dirac equation in spherical coordinates: 
Application of the multipole expansion of the electron— nuclei interaction 
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A non-perturbative approach to the solution ol the time-dependent, two-center Dirac equation 
is presented with a special emphasis on the proper treatment of the potential of the nuclei. In 
order to account for the full multipole expansion of this potential, we express eigenfunctions of the 
two-center Hamiltonian in terms of well-known solutions of the "monopole" problem that employs 
solely the spherically-symmetric part of the interaction. When combined with the coupled-channel 
method, such a wavefunction-expansion technique allows for an accurate description of the electron 
dynamics in the field of moving ions for a wide range of internuclear distances. To illustrate the 
applicability of the proposed approach, the probabilities of the K— as well as L- shell ionization of 
hydrogen-like ions in the course of nuclear a-decay and slow ion-ion collisions have been calculated. 

PACS numbers: 31.30.Jv, 34.80.Dp, 34.50.Fa 



I. INTRODUCTION 

Recent developments in accelerator and storage ring 
technologies have made it possible to perform a new 
generation of experiments on collisions between heavy, 
highly-charged ions. Of special interest in these studies 
are the low-energy collisions leading to the formation of 
short-lived quasi-molecular systems in which electrons 
move in the Coulomb field of two (or more) nuclei. Anal- 
ysis of the excitation, ionization, charge-transfer, and 
pair-production processes in such a low-energy domain 
may reveal important information about the properties 
and behaviour of few-electron systems and even of the 
quantum vacuum in the presence of extremely strong 
electromagnetic fields. To achieve and exploit the strong- 
field regime, a broad research program is planned to be 
undertaken at the future Facility for Antiproton and Ion 
Research (FAIR) in Darmstadt, at which ions up to bare 
uranium will be produced and decelerated to required 
energies [HQ. 

In order to better understand the basic atomic pro- 
cesses accompanying slow ion collisions, the experimental 
findings have to be supplemented by a detailed theoret- 
ical analysis. In the simplest case of the collision be- 
tween bare and hydrogen-like heavy ions, such an analy- 
sis can be traced back to the single-electron two-center 
Dirac problem. For small relative velocities and compa- 
rable charges of the nuclei, Z\ ~ Z%, the non-perturbative 
treatment of such a problem is usually required and can 
be performed by using various coupled-channel tech- 
niques. Along this line, the time-dependent electron 
wave packet is expanded in terms of eigensolutions of 
the stationary Dirac equation, which describes the two- 
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center system at a fixed internuclear distance R. The 
performance of the coupled-channel methods depends, 
therefore, on the efficiency of the spectrum generation of 
the time-independent Hamiltonian at each required R. 

An accurate solution of the static two-center problem 
is in general a rather sophisticated task which can benefit 
from a proper choice of coordinate system. During the 
last two decades in particular, a number of theoretical 
methods have been developed which make use of Cassini 
[H-IU and prolate spheroidal [1, 0] coordinate systems. 
Even though these (non-spherical) coordinates are very 
practical for the computation of quasi-molecular spectra 
at arbitrary internuclear distance R, their employment 
may be hampered by the lack of established numerical 
techniques for the evaluation of two-center matrix el- 
ements. Consequently, retention of standard spherical 
coordinates for the treatment of ion-ion (or ion-atom) 
collisions still attracts much current attention. The use 
of these, essentially one-center, coordinates for the de- 
scription of the two-center problem also requires the 
development of various approximate methods. Within 
the linear combination of atomic orbitals (LCAO) ap- 
proaches (8l4l2j. for example, quasi-molecular wavefunc- 
tions are constructed from sets of atomic orbitals, cen- 
tered on each nucleus. Yet another and very promising 
method relies on the direct solution of the two-center 
Dirac problem. Such a solution is rather straightfor- 
ward and well-elaborated if the electron-nuclei poten- 
tial is approximated by its spherically symmetrical part 
[l3l - fl7j ]. This so-called monopole approximation is suc- 
cessfully used for the description of strong-field phenom- 
ena in close-ion collisions, but performs poorly when the 
Coulomb centers are far from each other. The extension 
of the multipole theory towards accounting for higher 
terms in the decomposition of the two-center potential 
is crucial, therefore, for the proper treatment of heavy- 
ion collisions in spherical coordinates. 
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A number of efforts have been focused in the past on a 
straightforward solution of the (radial) Dirac equation for 
the complete two-center potential [3, [HI ■ In these stud- 
ies, the components of the quasi-molecular wavefunctions 
were found upon integration of an infinite system of cou- 
pled differential equations which account for all terms of 
multipolc expansion. An alternative and computation- 
ally very efficient approach to the two-center problem 
in spherical coordinates is proposed in the present work. 
We show that solutions of the stationary Dirac equation 
can be constructed for each internuclear distance R by 
means of a two-step procedure. As will be discussed 
in Section III Al the use of the dual kinetically balanced 
(DKB) B-spline basis set method [2(| for finding eigen- 
functions of the monopolc Hamiltonian constitutes the 
first step of the procedure. Based on the "monopolc" 
basis set, which is, thanks to the DKB algorithm, free 
of spurious, non-physical solutions, we generate then, in 
the second step, the two-center wavefunctions for any 
required number of multipoles in the potential expan- 
sion. Since the effective solution of the stationary two- 
center problem is, by itself, only an intermediate stage 
in the treatment of the time-dependent Dirac equation, 
the evaluation of the wave packet describing the electron 
dynamics in the field of moving nuclei will be discussed 
in Section III B[ In particular, we obtain the decomposi- 
tion of such a packet in terms of (stationary) two-center 
wavefunctions and determine the expansion coefficients. 
Although the developed approach can be applied to any 
collision between bare and hydrogen-like ions, indepen- 
dent of their charges and impact parameter, here we re- 
strict our analysis to two case studies of the electron 
loss in the course of (i) nuclear a decay and (ii) charge- 
symmetric ion-ion scattering at zero impact parameter. 
The first of these processes may be understood reason- 
ably well within the framework of the first-order pertur- 
bation theory [2ll - |2"H which will be employed in Sec. II V Al 
for testing the accuracy of our (non-perturbative) calcu- 
lations. In contrast, the ionization accompanying slow 
collisions between two heavy ions provides an example of 
a purely non-perturbative problem. To demonstrate the 
potential of the proposed method for tackling this prob- 
lem, we present in Sec. II V Bl predictions for the if-shell 
ionization in the U 91+ -U 92+ scattering. Based on the 
calculations conducted, we confirm a good performance 
of our time-dependent non-perturbative approach, pro- 
vided that the full multipolc expansion of the electron- 
nuclei interaction is taken into account. Summary of our 
results and a brief outlook will be given in Section [V] 
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II. THEORETICAL BACKGROUND 

The electron dynamics in the Coulomb field of two nu- 
clei is described by the time-dependent Dirac equation: 
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where the Hamiltonian reads, in spherical coordinates, 
as: 

H TC =a-p + V(Z u |r-Ri|) + V(Z 2 , |r-R 2 |) + /3. (2) 

In this expression, p = — zV is the electron momentum 
operator, (3 and a = {a x , a y ,a z } are the standard Dirac 
matrices, and the potential generated by the ith nucleus: 
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is a function of its charge density distribution p{r, Zj) 
and charge Z^. Moreover, Ri and R 2 describe positions 
of the nuclei with respect to the center-of-mass of the 
system: 
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where the internuclear vector R = R(t) varies over time. 

In what follows, we shall discuss the solution of the 
time-dependent Dirac equation (H])-© for relative ion 
velocities that are much smaller than the bound elec- 
tron velocity v w aZ,; . For such a slow collision regime, 
the adiabatic approach is justified and requires first the 
treatment of the static two-center problem. In the next 
subsection, therefore, we will show how the eigensolu- 
tions of the time-independent (two-center) Hamiltonian 
can be efficiently generated for any internuclear distance. 



A. Stationary two— center Dirac problem 

For each (instantaneous) position of the nuclei, the 
spectrum of the two-center system can be obtained by 
solving the time-independent Dirac equation: 



(5) 



where E is the total energy and the Hamiltonian Htc 
is given by Eq. ©. Analysis of such an eigcnproblcm, 
can be significantly simplified by the proper choice of the 
quantization (z-) axis. For example, by setting this axis 
along the internuclear vector 12, we can write the multi- 
pole expansion of the two-center potential from Eq. © 
in the form: 

F T c(r,R) = V(Z 1 ,\r-R 1 \)+V(Z 2 ,\r--R 2 \) 

OO 

= J2 V i( r > R ) P i( cos9 )> (6) 
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where Pi is the Lcgcndrc polynomial, 6 is the polar angle 
of the vector r, and the expansion coefficients V\ are given 
by: 



V t (r,R) 
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Moreover, if summation over I in Eq. ([7]) is restricted 
to the zeroth term, I = 0, the electron-nuclear interac- 
tion is governed by the spherically symmetric potential 
Vrc( r , R) = Vo ( r : -R)- The solution of the Dirac equation 
within such a monopole approximation is well-elaborated 
and has been discussed in a number of works [l3T[l7l | . In 
particular, the eigenfunctions of the monopole Hamilto- 



nian H, 
fori 



(0) 

TC 



a. ■ p + Vq (r, R) + /3 can be found in the 



'Hfl 



tf«(r)x- ¥ (f) 



(8) 



where Xkh is the standard Dirac spinor and the radial 
components satisfy the equation: 
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In order to solve this radial eigenproblem, we use the dual 
kinetically balanced (DKB) B-spline basis set method 
[20j . Since such a DKB approach has been widely applied 
in the past for the treatment of spherically symmetric 
Dirac problems, we will not discuss its details here. In- 
stead, we just mention that the DKB method allows one 
to avoid the spurious (non-physical) solutions of H^c 
and to generate a quasi-complete set of wavefunctions 
{</>™ M (r)}, n = 1, ...N for each value of the Dirac an- 
gular quantum number k. These functions describe the 
electron states with the energy e rlK in the spherically- 
symmetric potential Vo (r, R) , and their overall number 
N depends on the size of the basis set. 

The sets of eigenfunctions </>™ M (r), derived for the 
spherically symmetric problem can be employed to de- 
scribe the electron dynamics for relatively small distances 
between colliding nuclei. If R increases, the monopole ap- 
proximation is no longer valid and one has to account for 
the full two-center potential ([5]) when solving the time- 
independent Dirac equation ([5])- In the present work, we 
propose to present solutions of such an exact eigenprob- 
lem in terms of the monopole functions: 
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where K is a parameter limiting the number of partial 
waves in the sum. The expansion coefficients can be 



determined then based on the principle of least action, 
SS = 0, where the action is defined as: 
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By inserting the wavefunction (|10p into this expression 
and by evaluating the variation SS with respect to the 
change of expansion coefficients C*^ , we obtain a system 
of differential equations: 



OS 



0. 



which can be re-written in matrix form as follows: 
HC = EC. 

Upon evaluation of the elements of the matrix %: 

IK 



Ri,k = CkSi.k + 



J2Vi(r,R)P l (cos ( 



!=1 



OA 



(12) 



(13) 
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Eq. (|T3|) allows one to determine the vector C = 
{Ci, C2, •■••CV max }- Here, for the sake of brevity, we use 
short-hand notations C, 



^ and ej 



As seen from the discussion above, the spectrum of the 
time-independent Hamiltonian ([5]) for each fixed internu- 
clear distance R can be generated by means of the two- 
step procedure. In the first step, we employ the DKB 
finite-basis set approach to find solutions {e nK , </>™ M } of 
the monopole Hamiltonian. These solutions are used 
then, in the second step, to solve the generalized eigen- 
value problem (|13[) and to obtain both the expansion co- 
efficients of the wavefunctions $fc M (r) and the en- 
ergies of the electron states in the full two-center 
potential ©. In the next Section, such a new set of 
cigenstates {Ek^^kp,} will be employed for solving the 
non-stationary Dirac problem. 



B. Time— dependent two— center Dirac problem 

Having generated a (quasi-) complete set of cigen- 
states of the two-center Hamiltonian (|2|) at each in- 
ternuclear distance R, we are ready now to solve the 
time-dependent equation (TTJ) using the coupled channel 
method. Within this approach, the electron wavepacket 
^(r, t) is expanded: 



*(r,t) = 5^o fcM (t)*fc M (r,t) > 

kfj, 



(15) 



in terms of the functions &kp, which parametrically de- 
pend on the internuclear distance and, hence, on time t. 
In Eq. ([TBI , moreover, cifc M (£) are the time-dependent ex- 
pansion coefficients whose squares, ^^(t)) 2 , provide the 
occupation probabilities of the states |$fe M ) at a partic- 
ular instant in time. In order to find these coefficients, 
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we substitute the expansion (|15[) into the Dirac equation 
([1]) and derive the system of coupled channel equations: 



i-^j.a>kn{t) = E kl _ L {t)a kl _ L (t) 

-i 53 a nfl ,(t) | 



Any further analysis of this system requires the knowl- 
edge of how the electron-nuclei potential (J6j) varies with 
time. Since the time-dependence enters into the prob- 
lem solely through the intcrnuclcar distance R, the equa- 
tion of motion of colliding nuclei must be established. In 
the present work we consider the simplest case of motion 
along the Rutherford trajectories. In this case the time, 
the internuclear distance and the tilt angle of the molec- 
ular axis can be expressed in terms of the dimensionlcss 
parameter £ as follows: 



t 

R 



= — (esinh£ + £) 
= a(ecosh£ + 1) , 

= 2 arctan 



Ve2^T(tanh(C/2) + 1) 



e+1 - (e- l)tanh(£/2) 



(17) 
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with b denoting the impact parameter, Voo the asymp- 
totic value of the relative velocity of two particles at 
t = co, and M 12 the reduced mass. By inserting Eq. (fT7|) 
into the system of coupled channel equations (|16l) and 
re- writing it in terms of the parameter £, we derive: 
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where we used the relation 
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which is valid if the collision occurs in XZ-plane. Here 
j y is y-component of the total momentum projection op- 
erator and 8k n is the anti-Kronecker delta symbol. The 



paramctrization of t, R and 9 in terms of £ is most natural 
since the differential equation governing the time evolu- 
tion of R is autonomous i.e. an exact solution is possible 
only for t(R) and not, as required, for R(t). 

In order to solve the system of coupled channel equa- 
tions and, hence, to find the expansion coefficients a k , it 
is convenient to re-write Eq. (|19p in matrix form: 



(21) 



i^a(0 = M(£)a(£), 



where a = {ai,a 2 ,...}, and the individual elements of 
-Mfc,n(0 are given by 



M kfl! nil . (£) 
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The matrix equation (|21[) can be integrated numerically 
on a grid of spacing A£ according to: 



a(£ + A£) 



-iAf(M 



a(£) + 0(A£ 3 ), (23) 



and determines the vector a(£ + A£) at the "time" £ + A£ 
provided that the expansion coefficients coefficients a kfJl 
at the earlier moment £ are known. Since the matrix ex- 
ponential in the right -hand-side of Eq. (|23[) is unitary, 
the norm of the vector a will be preserved at each itera- 
tion. 

The iteration scheme (f2"3"|) represents the final step 
in the numerical treatment of the time-dependent two- 
center Dirac equation (JXJ) . In Section ITVl we will use this 
scheme in order to investigate the electron ionization in- 
duced by the nuclear a decay as well as the slow ion-ion 
collisions. In the present calculations, we shall restrict 
ourselves to the simplest case of zero-impact-parameter 
picture, 6 = 0. Within this framework, the last term 
in Eqs. (fT9f and (|20[) vanishes and, hence, the matrix 
elements of the evolution matrix (|22j) are diagonal in /x. 



III. DETAILS OF COMPUTATIONS 

Having discussed the non-perturbative approach to 
the solution of the two-center Dirac problem, we are 
ready now to investigate the electron emission accompa- 
nying both, the a-decay of heavy nuclei and the slow ion- 
ion collisions. Before starting with the presentation and 
analysis of the numerical results, let us briefly summarize 
the most important details of our calculations which, as 
mentioned in Section|TTl can be split into three stages. In 
the first step of this procedure, the eigenfunctions of the 
spherically-symmetric Hamiltonian H^c arc obtained by 
the DKB B-spline basis set method which guarantees the 
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absence of the non-physical spurious states in the spec- 
trum 20]. In the present work, we used about 200 B- 
splincs of eighth order defined in a box of size L ~ 10 5 fm 
in order to construct "monopole" wavefunctions 0™ M (r) 
with energies in the range < e nK < 10 mc 2 . Based 
on the detailed numerical analysis, we argue that such a 
truncated basis set allows one to achieve ~ 5-10 % ac- 
curacy in the prediction of the ionization cross sections. 
As the second step of the non-pcrturbative treatment, 
the solutions of the full two-center Hamiltonian J2]) are 
expanded in terms of </>™ M (r) (cf. Eq. I|10p ). Along this 
line, we obtain about 300 functions $fc M (r;i?) and cor- 
responding energies Ek for each internuclear distance R 
(or, equivalently, dimensionless parameter £ (|17|) ). It is 
worth mentioning that the solutions of cigcnproblcm (|13[) 
and, hence, $fc p (r) are defined up to an arbitrary sign. 
In our calculations, this sign is chosen for all ^^(r) from 
the requirement that their large radial components, cal- 
culated for two successive steps over £, behave similarly 
near the origin of the coordinates, i.e. for r = 0...500 fm. 

With the help of generated basis sets 
{^k/j,{ r 'i R(Q)}k=i,...N we are finally able to per- 
form the time propagation of the electron wavepacket in 
the field of moving nuclei. Prior to starting this prop- 
agation, one has to define the electron wavefunction in 
the initial moment of time. Indeed, the initial conditions 
depend on the particular process under consideration. 
For the nuclear a-decay, for example, we assume that 
the electron is originally in the ground Is 1/2 state of the 
united nucleus of charge Z. Since the time-propagation 
begins from the moment when the a particle leaves the 
potential barrier at the distance Ro ~ 10 fm from the 
daughter nucleus, we project, at £ = (corresponding 
to t = 0), the wavefunction ipi Sl/2 (r;Z) onto the basis 
set of eigenfunctions of Hamiltonian ^ describing 
system of two Coulomb centers with charges Z\ = 2 
and Z2 = Z — 2, placed at distances R\ and R2 with 
respect to their ccnter-of-mass (see Eq. Such 
a projection procedure allows us to account for the 
shake-off effect and to obtain the first set of expansion 
parameters {<Xfc M (£ = 0)}fc = i r ..jv which are used then to 
find the electron wavepacket in subsequent time steps 
(cf. Eq. (|2"3")l ). In order to produce results, presented 
in the Section IIV Al time propagation was carried out 
for about 750 such steps of A£ = 0.01, this corresponds 
to the retreat of the a particle to a distance of about 
10 4 fm. 

In contrast to the a-decay, the time propagation of 
the electron wavepacket in the field of two colliding ura- 
nium ions, studied in Section HVBl was started from the 
moment when the ions are separated from each other by 
the distance R = 5 ■ 10 3 fm. In this initial moment, 
the electron finds itself in the ground lsi/2 state of one 
of the projectiles. The wavefunction of such a state is 
given by the sum of the lowest-lying gerade and unger- 
ade solutions of the (stationary) two-center Hamiltonian, 
ipu 1/a (ri; Z = 92) w l/-v/2($i«7 S + $ i(Ju ); anapproxima- 
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FIG. 1: Tf-shcll ionization probability of hydrogen-like xenon 
(top panel), gadolinium (middle panel) and polonium (bottom 
panel) ions following the a-decay. Non-perturbative calcula- 
tions were carried out within the monopole approximation 
(dashed line) and by taking the full two-center potential into 
account (solid line). The probability is scaled xlO . 



tion whose quality increases with the number of partial 
waves in the expansion ([TU)) . In the calculations bellow 
all partial waves with the Dirac angular quantum num- 
ber in the range k = — 10... + 10 arc employed leading to 
about 10 % accuracy of the presented results. 



IV. RESULTS AND DISCUSSION 

A. Ionization following a— decay of heavy nuclei 

The non-pcrturbative approach presented in Sectionim 
can be used to study basic atomic processes accompany- 
ing slow collisions of two ions independent of their nu- 
clear charges Z\ and Zi- In this section, we employ it 
to re-analyze the nuclear a decay, which is an example 
of a (charge-) asymmetric collision, Z\ « Z^, with zero 
impact parameter and which can be treated also within 
first-order perturbation theory. As mentioned already, 
such a perturbativc treatment has been successfully ap- 
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plied over the last decades in a large number of studies 
j22l - |24l . [27j . In order to compare predictions of the non- 
perturbative and perturbative theories, we consider the 
decay of a-active 110 Xe, 148 Gd and 210 Po isotopes. For 
these zero-nuclear-spin nuclei, calculations have been 
performed for the ionization probability Pk of an electron 
from the ground lsi/2 state of an initially hydrogen-like 
system. In Fig. [TJ we display the non-perturbative results 
for the Pk as a function of the internuclear distance. To 
deduce this probability, we have evaluated the electron 
wavefunction ^(r, t) at each step of the time propagation 
(see Eqs. (fl"5)) - ([2"3"| and related discussion) and projected 
it onto the positive-energy solutions of the two-center 
Dirac equation for the instantaneous distance R = R{t) 

P K (R(t)) = Yl l(^(r)|*(r,t))| 2 

E k >mc 2 

= E \ a *M 2 - (24) 

E k >mc 2 



TABLE I: K-sheYl ionization probability of hydrogen-like 
xenon, gadolinium and polonium ions following the a-decay. 
The non-perturbative calculations, performed for R — > oc 
by using the monopole as well as exact approximations to 
the two-center potential, are compared with the first-order 
perturbation results a nd p redictions by Law [2^ ] . and Fis- 
chbeck and Freedman [2y]. The asymptotic kinetic energy 
of a particle Tki n = from Ref. [2^| is given in the 

second column. All probabilities are of the order of xlO 6 . 



Ion 


Tkin 


Perturbative 


Non-perturbative 




(MeV) 




monopole exact 


HU Xe +b3 


3.7 


3.61 


2.6 3.2 


148 Gd+ 63 


3.1 


2.15 


1.6 2.3 


210p Q +83 


5.4 


2.00 


1.4 2.1 






1.81" 








2.03 6 





"Law [32] 

''Fischbcck and Freedman [2^| 



Calculations have been performed both within the 
monopole approximation, in which summation over I in 
Eq. (JB|) is restricted to the zeroth term, and by taking the 
full two-center potential Vrc( r , R) into account. As was 
expected, these two approaches agree only for relatively 
small internuclear distances. If R becomes greater than 
500 fm, the monopole approximation can significantly 
underestimate the ionization probability; an effect which 
becomes most pronounced for the heavy nuclei. 

Fig.[T]shows that at very large distances, R > 8000 fm, 
the ionization probability Pk converges to some final 
value which depends only on the charge of the mother 
nucleus and the initial velocity of the a particle. This 
"asymptotic" value of Pk is displayed in Table Q] for 
xenon, gadolinium and polonium ions, and compared 
with the results of our first-order perturbation calcula- 
tions (see Ref. (24) for further details). Moreover, the 



previous (perturbative) predictions of Law [221 ] , and Fis- 
chbcck and Freedman [29| obtained for the decay of polo- 
nium are given in the third column. As seen from the 
table, the non-perturbative treatment, based on the full 
multipole expansion of the two-center potential, repro- 
duces well the ionization probabilities for all three ions. 
In particular, both perturbative and non-perturbative 
theories yield results that agree to within 5 % if applied to 
the exploration of the a decay of polonium ions. If, how- 
ever, the potential Vrc( r 5 R) is approximated in Eq. @ 
by the single monopole term, the non-perturbative calcu- 
lations may result in approximately a 30 % misestimation 
of Fx- 

Until now we have discussed the a-decay-induced ion- 
ization of hydrogen-like ions that have been prepared ini- 
tially in the ground lsi/2 state. In order to verify the 
performance of the non-perturbative technique, based on 
the multipole expansion of the two-center interaction op- 
erator, it is also worth considering the electron emission 
from the various L subshclls. Even though experimen- 



tal observation of the L-shell ionization of hydrogen- 
like systems might be hampered by the short lifetimes 
of excited ionic states, it can be measured for neutral 
atoms. Theoretically, such an atomic inner-shell ion- 
ization can be well described by using the developed 
approach if the proper screening potential is used in 
Eq. <j3j) . The analysis of the screening effects in a-decay- 
induced processes in neutral systems is, however, out of 
the scope of the present work. Instead, we just employ 
the L-shcll ionization of hydrogen like ions as a test- 
ing ground for the non-perturbative theory from Sec- 
tion |lll The internuclear-distance-dependent probabili- 
ties for the ionization of 2si/2 (top panel), 2pi/ 2 (middle 
panel) and 2p 3 / 2 (bottom panel) states of hydrogen-like 
polonium are evaluated based on this theory and are pre- 
sented in Fig. [5J Similar to before, calculations have 
been performed by accounting for the full multipole ex- 
pansion of the two-center potential (solid line) and by 
restricting this summation to the monopole term only 
(dashed line). Agreement between these two approaches 
can be observed again only for small internuclear dis- 
tances, while for R > 600 fm the monopole calculations 
underestimate the ionization probabilities by more than 
25%. Moreover, the monopole approximation fails to re- 
produce i"2p 3 /2 f° r the entire range of R. 



The asymptotic values of P2s 1/2 , P2 Pl/2 aim r- 2p3/2 
eulated for large distances R are presented in Table lU 
and compared with the predictions of first-order pertur- 
bation theory [24J and data by Law [22j ■ As in the case of 
JT-shell ionization, the full account of the electron-nuclei 
interaction Vrc{ r i'R-) in Eq. (0) leads here to approx- 
imately 5% agreement between the predictions of per- 
turbative and non-perturbative theories for the entire 
L shell. In contrast, the time propagation of the elec- 
tron wavepacket in the spherically-symmetric potential 
Vo(r, R) yields the probabilities Pl that are 30% smaller 



and P?. 



cal- 
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FIG. 2: Ionization probability of the 2s\/2 (top panel), 
2pi/2 (middle panel), and 2j>3/ 2 (bottom panel) states of 
hydrogen-like polonium Po following the a-decay. 
Non-perturbative calculations were carried out within the 
monopole approximation (dashed line) and by taking the full 
two-center potential into account (solid line). The probability 
is scaled xlO 5 . 



comparing to the perturbative results. Again, these find- 
ings stress the importance of the higher multipole contri- 
butions to the electron-nuclei interaction for the time- 
dependent analysis (|23p of the electron dynamics accom- 
panying ion collisions. 



B. Ionization in U 91+ -U 92+ collisions 

So far, we have shown that the time-dependent 
method (JTHJ), based on the expansion of the basis wave- 
functions in terms of monopole solutions, can be success- 
fully utilized to study the a-decay-induced ionization. 
Besides this — purely perturbative — problem, the perfor- 
mance of the developed approach has been also examined 
for slow collisions between two high-Z ions. In contrast 
to the a-decay, theoretical analysis of such collisions usu- 
ally can not be carried out within the framework of the 
perturbation theory and demands the application of non- 



TABLE II: L-subshell ionization probabilities of hydrogen- 
like polonium 210 p o + 83 following the a-decay. The non- 
perturbative calculations, performed for R — > oo by using the 
monopole as well as exact approximations to the two-center 
potential, are compared with the first-order perturbation re- 
sults and predictions by Law [z| . Probabilities given in units 
xlO and for the kinetic energy of the emerged a-particle 
T kin = 5.4 MeV 



State 


Perturbative 


Non-perturbative 
monopole exact 


2si/ 2 


4.80 


4.0 4.96 


4.75" 




2P1/2 


0.54 


0.30 0.64 


0.50 a 




2P3/2 


0.61 


0.04 0.61 


0.60° 





"Law [12 



perturbative techniques. Along this line we have focused, 
in particular, on the if -shell ionization in u 91 +-U 92 + col- 
lisions at zero impact parameter. The ionization proba- 
bility Pk has been calculated based on Eq. (|24|) . where 
the electron wavepacket ^(r, t) was propagated from a 
time when the ions were at a distance R = 5 ■ 10 3 fm, 
through the closest approach R$ w 50 fm, to a mo- 
ment when the internuclear distance increased again to 
R = 5 • 10 3 fm. In Fig. for example, Pk is displayed as 
a function of the distance R and for the (relative) colli- 
sion energies T p = 1.8, 2.0 and 2.2 MeV/u. As seen from 
the figure, the steep rise of the ionization probability ap- 
pears immediately after the point of closest approach R 
at which the (relative) ionic motion is suddenly reversed 
and the electron can be "shaken off" into the contin- 
uum. Such a behaviour of the Pk as well as its further 
damped oscillations have been predicted previously in 
Ref. [lH based on the monopole approximation and now 
is confirmed by our theory that accounts for the multipole 
expansion of the electron-nuclei interaction. Moreover, 
our calculations clearly indicate a rise of the ionization 
probability with the collision energy. For example, the 
asymptotic value of the Pk is increased by almost factor 
of three if the initial (relative) energy changes from 1.8 
to 2.2 MeV/u. Further significant enhancement of the 
Pk is predicted for higher energies at which the "div- 
ing" of the ground quasi-molecular state into the Dirac's 
negative continuum takes place (i~3j . However, since the 
analysis of such strong-field phenomena is out of scope 
of the present paper, we restrict here our calculations to 
the "undercritical" energy range, T p < 2.3 MeV/u. 



V. SUMMARY AND OUTLOOK 

In summary, we have laid out a theoretical approach to 
the time-dependent two-center Dirac problem. Within 
such an approach, the wave functions, describing the (sin- 
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FIG. 3: Ionization probability of the lsi/2 state of hydrogen- 
like uranium colliding with bare U 92+ ion. Non-perturbative 
calculations, based on the multipole expansion of the 
electron-nuclear interaction operator, were carried out for 
zero impact parameter and for the initial relative kinetic en- 
ergy of ions T p = 1.8 MeV/u (solid line), 2.0 MeV/u (dashed 
line), and 2.2 MeV/u (dash-dotted line). The negative and 
positive values of R correspond to the times when ions ap- 
proach and move away from each other, respectively. 

gle) electron dynamics in the field of two moving nuclei, 
are expanded in terms of solutions of the stationary Dirac 
equation (0. We have argued that these stationary so- 
lutions can be efficiently constructed in spherical coor- 
dinates and for each internuclear distance R by means 
of the two-step procedure. The first step of the pro- 
cedure consists of employing the dual kinetically bal- 
anced (DKB) B-spline basis set method to find eigen- 
functions of the Hamiltonian Hq , which accounts for the 
spherically-symmetric part of the electron-nuclei inter- 
action. On the basis of these functions we generate, in 
the second step, the required solutions of the stationary 
two-center problem. 

The developed time-dependent approach can help to 
explore various atomic processes accompanying slow ion 
collisions. In the present work, for example, we used this 
theory to calculate the electron-loss probabilities for the 
(i) a-decay of hydrogen-like xenon, gadolinium and polo- 
nium ions, and (ii) U 91+ -U 92+ scattering at zero impact 
parameter. The a-decay, being an example of charge- 
asymmetric collisions, can be described sufficiently well 
within the framework of first-order perturbation theory. 
Calculations based on this theory have been used to prove 



the accuracy of our non-perturbative approach. For the 
K- and L-shell a-decay-induced ionization, predictions 
of both perturbative and non-perturbative methods were 
found to agree to within about 5 % if the multipole ex- 
pansion of the two-center potential is taken into account 
in the time-dependent Hamiltonian @. If, in contrast, 
this potential is approximated by its monopolc term, our 
calculations may underestimate the ionization probabil- 
ities by more than 30 %; this failure of the monopolc 
approximation becomes most pronounced for large inter- 
nuclear distances. Based on these findings we stressed the 
vital importance of the proper treatment of the electron- 
nuclei interaction for the accurate description of slow ion- 
ion collisions. The rigorous "multipole" approach has 
been employed then to explore the if -shell ionization 
accompanying TJ 91+ -U 92+ collisions. For this — purely 
non-perturbative — process, wc qualitatively confirmed 
the impact-parameter-behaviour of the ionization prob- 
ability, which was predicted previously by Betz and co- 
authors [l3| within the monopole theory. 

Both the a-decay of hydrogen-like heavy ions and the 
tj9i+_tj92+ sca ttering have been explored in the present 
work for the case of zero impact parameter. Of course, 
the developed non-perturbative method is not limited 
to such a simple geometry and can be applied to analyze 
heavy-ion collisions at b ^ 0. For these collisions, the last 
term of Eqs. (fT9)) - (j20"|) . that accounts for the rotation of 
the internuclear distance, does not vanish and makes the 
elements of the evolution matrix (|22p non-diagonal in /i. 
The impact-parameter dependence of the electron loss as 
well as the excitation and the charge transfer processes 
will be discussed in a forthcoming publication and will 
help in planning future experiments on slow collisions 
between two high-Z projectiles. These experiments are 
likely to be carried out at the Facility for Antiproton and 
Ion Research (FAIR) in Darmstadt and are expected to 
reveal unique information about the quantum electrody- 
namics of extremely strong fields. 
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